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• Premise of the study: Fifty micro satellite loci were developed for the endangered species Juglans mandshurica to investigate 
its genetic diversity and population structure. 

• Methods and Results: In all, 50 micro satellite markers were isolated from /. mandshurica, using the Fast Isolation by AFLP of 
Sequences Containing repeats (FIASCO) protocol. Twenty of these polymorphic markers were assessed in samples collected 
from 98 individuals among five populations in northeastern China. Across all of the /. mandshurica samples, the number of 
alleles per locus ranged from one to 17. 

• Conclusions: These new microsatellite loci will be useful for conservation genetics studies of /. mandshurica. 
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Juglans mandshurica Maxim. (Juglandaceae) is a deciduous 
tree that is widely distributed in northeastern China, as well as in 
some areas of Korea and the Russian Far East (Lu, 1982). How- 
ever, the habitat of /. mandshurica is now seriously threatened 
by human activities (Wang et al., 2011). Juglans mandshurica 
has been used in folk medicine for many years, particularly in 
China and India (Xu et al., 2010). Its leaves, fruits, roots, stem 
bark, and seeds have also been used as traditional medicine for 
cancer treatment in Asia and Europe (Kim et al., 1998; Sun, 
2004). The therapeutic benefits of /. mandshurica have been as- 
cribed to its naphthoquinone content. Naphthoquinone is known 
to have various physiological properties including the induction 
of apoptosis effect (Kang et al., 2001). 

To date, there have been only a few reports on microsatellites 
in /. mandshurica (Qi et al., 201 1). The molecular genetic diver- 
sity studies of /. mandshurica are limited to the work by Woeste 
et al. (2002) using dominant markers. Here, we used an im- 
proved technique for the isolation of codominant compound mi- 
crosatellite markers (Lian et al., 2006) to isolate 20 polymorphic 
microsatellite loci from /. mandshurica. These polymorphic 
markers can be useful for the development of effective conser- 
vation programs of the species. 
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METHODS AND RESULTS 

Total genomic DNA was extracted from silica gel-dried leaves of J. mandshu- 
rica using a modified cetyltrimethylammonium bromide (CTAB) method (Doyle, 
1991). After digestion with EcoRl and Hindlll (TaKaRa Biotechnology Co., 
Dalian, Liaoning, China) at 37°C for 4 h, a fraction containing 200-1000-bp frag- 
ments was isolated from total genomic DNA (250 ng) and ligated with EcoKL 
adapters (5'-CTCGTAGACTGCGTACC-3' and S'-CTGACGCATGGTTAA-S') 
and Hindlll adapters (5'-GACGATGAGTCCTGAG-3' and 3'-TACTCAGGA- 
CTCTCGA-5') using T4 DNA ligase (Biomanbio, Shanghai, China) at 4°C over- 
night. Ligated fragments were subsequently PCR amplified with adapter-specific 
primers. To enrich the fragments containing microsatellite repeats, PCR products 
were hybridized with biotin-labeled probes, (AG) 15 and (GT) 15 . Streptavidin-coated 
magnetic beads (Promega Biotech, Beijing, China) were then used to separate and 
capture DNA fragments that were hybridized to the probes, according to the manu- 
facturer's instructions. DNA fragments containing simple sequence repeats were 
eluted from streptavidin-coated magnetic beads in a 1.5-mL tube and put in a 
freeze-drying machine at -40°C for 3 h. Microsatellite-enriched fragments in 
eluted solutions were PCR amplified in 20- uL reaction mixtures, consisting of 
8.6 uL of DNA template, 10 uL of 2x PCR Mix (Boyouxinchuang Biotech, 
Beijing, China), 0.3 uL of E00 (10 uM), 0.3 uL of H00 (10 uM), and 0.2 uL of 
Taq polymerase. PCR conditions were as follows: five cycles of denaturation at 
94°C for 45 s, annealing at 50°C for 45 s, and extension at 72°C for 45 s; and 
final extension at 72°C for 10 min. The PCR products were ligated into the 
pMD-18T vector (TaKaRa Biotechnology Co.) and then transformed into 
Escherichia coli DH5a competent cells to generate the microsatellite sequence- 
enriched library. Recombinant colonies were selected by white/blue screening 
on Luria-Bertani agar plates containing ampicillin (60 ug/mL). 

Positive clones were then screened using the universal primers of pMD-18T 
(M13-47 and RV-M, 5'-CGCCAGGGTTTTCCCAGTCACGAC-3' and 5'-GAG- 
CGGATAACAATTTCACACAGG-3', respectively), as well as the (AG) 9 / 
(GT) 9 tandem repeat primers. PCR amplifications were performed in a 15-uL 
reaction mixture containing DNA template (single colonies), 7.5 uL of 2x PCR 
Mix (Boyouxinchuang Biotech), 0.1 uL of M13-47 (0.2 uM), 0.1 uL of RV-M 
(0.2 uM), 0.1 uL of (AG) 9 /(GT) 9 (0.2 uM), 0.2 uL of Taq polymerase, and 7 uL 
of ddH 2 0. PCR amplification was performed in the GeneAmp PCR System 9600 
Thermal Cycler (Applied Biosystems, Foster City, California, USA) under the fol- 
lowing conditions: initial denaturation at 94°C for 10 min; 33 cycles of denaturation 
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Table 1. Primer sequences and characteristics of 50 micro satellite loci in Juglans mandshurica. 



Locus 



Primer sequences (5'-3') 



Repeat motif 



Fragment size (bp) 



r a (°C) 



GenBank accession no. 



CCAAGGGAAATACAAGGTCT 

GTTGCCTGAACATCACAGAT 

GGTCAAGGCTCTCTGCCTCAA 

TAAAATCACCCCCCTCCACTC 

GGAGTTTCGGGTAGGGTTTGA 

TAATGGTTTGGAGGAATGGAG 

TAGGAGATATTTTCAAGAGGA 

TTGTGTGAGCATGAGTTTTAG 

TATTTCATTGCCAAGACCGAG 

GTCCCTAAAGAGTGATTTTGTTGT 

AAAC AGC AT C T GAAAC C AC A 

GCTTTTGCCTCCATAATTAG 

AAAAGAGGTGTTGAGGATGG 

TTATTTTCACTTGCTTTGCC 

C T GAAAAGT GGGGC AAAGC A 

GGGAGACATACCCGACAAGG 

CACGACGGCACAACTAAAGG 

CACTGAAGGCACACCCAAGA 

AAGATTGTTTCTAAGTTTGTC 

AATGTGTAGGTCAATAGAGG 

TCCACCGTAAAAGATTGTT 

GCTTCATAGAGATTCCCAT 

AAACCCTATCTTCCGCAGAA 

GATGGAGAGCTAAGGAGTCG 

GGTTAGAGTGAGCGAGAGTTG 

ATCCTTAGAGTTGAATGGGC 

AAAAAAC C T T GC G AC CAAAG 

C C AC T AAAAC AC AC T C C AT C A 

CTCTCTCGCAAAACCTCTCG 

CGTGCATGCTAGGAAACTTA 

CTGCTGATGTGGTGGTGAAG 

GATCTGGTTTGTGGAGGGAA 

TGAGTACAGAACTGGCATG 

ACTGGGGTGGGCTAAAAAG 

TAGAAACCTCGTGGACTTG 

ACTGCTAGAGCCTATGGAA 

CCGAATGGAGGAAGGAAGG 

ACACGGCTTAGGGCCATAAA 

CAGCCACCCCATTACCATTC 

C AC AAC T C AC AAAAC C AAC AAC 

TGTTACTCTGTTGGGTCGT 

CTGGTGTAGCAGTTCATTT 

CCTTCTGCT C AGAC AAAC A 

CTAGACCAAAGACCACAACTAT 

CAGGACACGACAACCCCAT 

TAACACTCCACGCACGCAC 

TATATTGTGGGAGGTGGGT 

TTAAGGGAGTTGTTGAAGC 

GGGAGAATGAAAATGACGG 

GAAAC GAAGGAAAAATGAGG 

CTTCAAATAATGGAACGGT 

TAAAGAGATGGAGTACGCT 

TTCATAGCACATACCAGTTC 

TCCGTAACATCAATCATTC 

GCTTTTGTTACTTTGTGCC 

GAGCTGAATTTTTTACCTGA 

GAAT AAAAGAAGT T T GAC 

TCTGACCCAAATTCCATAG 

ATGGGAATCACAGGTGAC 

TGGGAATATCTTCGCAGAG 

CCTTCAACCACTTCAATA 

GAGACACACGCACAAAACC 

AAGTAAAACCCTAAGTCC 

AAAATCCAAACTCAAGCCC 

CCAAGGCACCAACACAAT 

CCCATGCAATAACAAACCA 

TTCCCCATCTCAAATCTG 

TTCCTCGTGTAAGTACCCC 

TGCTTTTACCTTCCCTCT 

GTCCAACCCAAGTTCCTCC 



(AG)„ 

(TQ 15 

(TQ 15 

(AC) 8 
(CA) 8 (AT) 5 (AG) 
(CT) 10 (CA) 8 T 

(GA) 13 

(GA) 14 

(AG) 20 
(AT) 3 (AG) 15 

(AG) 16 

(CT) U 

(AG) 15 

(GA) 26 

(CT) 26 

(AG) 9 

(TG) 18 

(AC) 10 

(GT) 22 
(CT) 9 (CA) 7 
(GA) n (TA) 4 

(AG) 13 
(AG) 3 (AG) 10 

(CT) 9 

(AG) 19 

(CT) 17 

(CT) n 

(CT) U 

(CT) 10 

(AG) 13 

(CT)i3 

(CT) 14 

(TC) 15 

(GA) 17 

(AG) 18 



186 
249 
184 
136 
164 
288 
113 
280 
144 
136 
108 
108 
186 
211 
171 
194 
187 
214 
126 
283 
196 
241 
200 
316 
124 
195 
283 
196 
191 
252 
244 
241 
266 
159 
296 



58 
57 
58 
57 
57 
56 
55 
59 
57 
58 
56 
57 
58 
56 
58 
59 
55 
58 
56 
57 
58 
56 
59 
57 
57 
56 
58 
57 
58 
57 
55 
57 
56 
58 
59 



JQ618127 
JQ618128 
JQ618129 
JQ618130 
JQ618131 
JQ618132 
JQ618133 
JQ618133 
JQ618134 
JQ618135 
JQ618136 
JQ618137 
JQ618138 
JQ618139 
JQ618140 
JQ618141 
JQ618142 
JQ618143 
JQ618144 
JQ618145 
JQ618146 
JQ618147 
JQ618148 
JQ618149 
JQ618150 
JQ618151 
JQ618152 
JQ618153 
JQ618154 
JQ618155 
JQ618156 
JQ618157 
JQ618158 
JQ618159 
JQ618160 
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Table 1. Continued. 



Locus 




Primer sequences (5'-3') 


Repeat motif 


Fragment size (bp) 


TA°Q 


GenBank accession no. 


P36 


F 

R 


AC C AAAC AG AAC G AG T AA 
AGGGGATAGATTGTGATAC 


(GA) 13 (GT) 4 


240 


57 


JQ618161 


P37 


F 

R 


CGTTTTGAGTTTCTGCCT 
CACACACTACATGGATGTC 


(GA) 33 


291 


58 


JQ618162 


P38 


F 

R 


ATGTGTAACCAGATAAGG 
ACAAATTGGTCATCCTAGA 


(GA) 18 


158 


57 


JQ618163 


P39 


F 

R 


GATCTGAAGAGCCTGCCT 
ATCCACCAAAACCCTAAAA 


(AG) 30 


198 


56 


JQ618164 


P40 


F 

R 


CCCTCGTTACTCTCCCTT 
ATGTGGGTCGTGGGTTTGTC 


(CT) 24 


170 


58 


JQ618165 


P41* 


F 

R 


AGCATACCTTCAATGGAT 
CGCAGAATACACGCAAATAG 


(AG) 18 


206 


59 


JQ618166 


P42 


F 

R 


ACAGAATGTGAGGTTTACTACG 
GTGCTACTCTTTTTGTGGAT 


(CT) 10 


250 


57 


JQ618167 


P43 


F 

R 


TACAGCACACCCCTGAAAT 
AACGACGCCGACCAACAAC 


(AG) 18 


234 


58 


JQ618168 


P44* 


F 
R 


AAGCATCATCTCTATTTCTC 
ACTTTGTGGGTGTTTCTAT 


(CT) n 


200 


57 


JQ618169 


P45* 


F 
R 


TTGATGCCTGTAGGTGAATG 
TAATGCTATGGAAGTATGGA 


(CT) 10 (CA) 13 


252 


57 


JQ618170 


P46 


F 

R 


ATTCTAATGCCGTCACTTG 
AGGATCTGTTGACACATACAT 


(ACC) 3 (AC) 8 


234 


58 


JQ618171 


P47* 


F 
R 


TGT G AAAAGG AC T T C AC AT 
GGATCTCAAGACTGGCTAG 


(AT) 4 (AG) 2 (AC) 


217 


57 


JQ618172 


P48 


F 
R 


TTATTATGATGGGGTCTTTG 
TTGATACCAACTGTAACGC 


(TG) 20 


198 


56 


JQ618173 


P49 


F 
R 


CAAAGCAGAGGCTGAGATA 
CTAAGCAAAGGAACGAAAG 


(GC) 4 (AC) 5 


180 


57 


JQ618174 


P50 


F 
R 


GGGGAAACAGCAACAATAG 
TCCAGGAAGTTAGGGTGAG 


(AC) 18 


268 


58 


JQ618175 



Note: T a = annealing temperature when run individually. 

* Indicates polymorphic primers; all other loci were monomorphic. 



at 94°C for 45 s, annealing at 58°C for 45 s, and extension at 72°C for 45 s; and 
final extension at 72°C for 10 min. A total of 90 colonies were randomly selected 
and sequenced on an ABI Prism ABI 3730xL automated DNA sequencer 
(Applied Biosy stems); 82 of these contained repeats. Primer Premier 5 software 
(PREMIER Biosoft International, Palo Alto, California, USA) was used to 
design primer pairs for 50 of the 82 repeat-containing samples. The forward 
polymorphic primer of each pair was also labeled with a fluorochrome, 6-FAM 
(Applied Biosystems). 

These 50 primer pairs were then tested for amplification in 20 samples from 
five different populations (four individuals each) from the mountains of north- 
eastern China (Appendix 1). PCR products were electrophoresed on 8% dena- 
turing poly aery lamide gels and visualized by silver staining. We found that all 
50 primer pairs successfully amplified products after PCR optimization (dena- 
turation at 94°C for 45 s, at 55-60°C for 45 s, and 32 cycles at 72°C for 45 s, 
with a final extension of 10 min at 72°C) (Table 1), and that 20 of these primer 
pairs amplified high-quality, polymorphic PCR products. These 20 polymor- 
phic primer pairs were characterized in 98 individuals among five natural popu- 
lations of J. mandshurica (Appendix 1). The distances between these individuals 
were more than 100 m. PCR products were analyzed by capillary electrophoresis, 
visualized using the ABI 3730xL sequencer with a GeneScan 500 ROX Size 
Standard (Applied Biosystems), and scored using GeneMarker (SoftGenetics, 
State College, Pennsylvania, USA). 

Genetic diversity statistics, including the number of alleles, observed 
heterozygosity, and expected heterozygosity, were calculated based on 98 sam- 
ples using the software POPGENE32 (version 1.31; Yeh et al., 1999). The 
number of alleles ranged from one to 17 (mean = 6), whereas the expected and 
observed levels of heterozygosity were 0-0.925 and 0-1.000 (mean = 0.358 
and 0.627), respectively. The large number of alleles per locus indicates the 
potential usefulness of these primers to characterize the population genetic 
structure of J. mandshurica. Moreover, nine loci (P8, P12, PI 8, P23, P26, P30, 
P41, P44, P45) in J. mandshurica were found to show significant deviations 
from Hardy-Weinberg equilibrium (P < 0.01 in % 2 test) (Table 2) due to 
heterozygote deficiency. No significant linkage disequilibrium was detected 
between any pairs of loci. 



CONCLUSIONS 

Twenty polymorphic microsatellite markers were identi- 
fied in /. mandshurica and used to investigate its genetic di- 
versity and population structure. The microsatellite loci 
described here (including both monomorphic and polymor- 
phic primers) can be useful for conservation genetic studies 
of J. mandshurica. 
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Table 2. Results of PCR screening of 20 polymorphic micro satellite loci in Juglans mandshurica.* 



Xiaoxing'anling Zhangguangcailing Changbai Mountain Laoyeling Wanda Mountain 



Locus 




(N= 19) 






(#=20) 






(#=20) 






(N= 19) 






(#=20) 




HWE P value b 


A 


H 0 


H e 


A 


H 0 


H e 


A 


H 0 


H e 


A 


H Q 


H e 


A 


H 0 


H e 


PI 


5 


0.313 


0.720 


6 


0.688 


0.470 


6 


0.250 


0.728 


5 


0.250 


0.754 


6 


0.125 


0.758 


0.128 


P3 


5 


0.250 


0.780 


5 


0.625 


0.474 


1 


0.375 


0.798 


5 


0.438 


0.760 


6 


0.250 


0.841 


0.012 


P5 


5 


0.375 


0.726 


6 


0.000 


0.712 


8 


0.125 


0.758 


5 


0.188 


0.601 


6 


0.125 


0.677 


0.011 


P8 


1 


0.438 


0.798 


4 


0.188 


0.591 


8 


0.250 


0.768 


6 


0.063 


0.750 


7 


0.313 


0.802 


0.000 


P12 


5 


0.500 


0.595 


3 


0.063 


0.619 


7 


0.563 


0.633 


6 


0.313 


0.534 


3 


0.625 


0.542 


0.000 


P14 


12 


0.400 


0.869 


8 


0.267 


0.770 


9 


0.438 


0.605 


8 


0.125 


0.794 


9 


0.077 


0.846 


0.010 


P16 


5 


0.188 


0.565 


3 


0.063 


0.542 


5 


0.500 


0.559 


4 


0.438 


0.569 


4 


0.375 


0.595 


0.113 


P18 


5 


0.688 


0.750 


4 


0.688 


0.466 


7 


0.750 


0.651 


5 


0.438 


0.655 


7 


0.375 


0.813 


0.000 


P21 


7 


0.125 


0.802 


6 


0.500 


0.629 


11 


0.125 


0.865 


8 


0.188 


0.728 


7 


0.250 


0.710 


0.011 


P23 


3 


0.563 


0.365 


4 


0.063 


0.625 


2 


0.625 


0.315 


3 


0.688 


0.365 


2 


0.438 


0.466 


0.000 


P26 


8 


0.188 


0.843 


5 


0.125 


0.649 


6 


0.063 


0.790 


5 


0.000 


0.665 


7 


0.000 


0.702 


0.000 


P28 


7 


0.313 


0.633 


4 


0.813 


0.236 


6 


0.125 


0.748 


4 


0.625 


0.375 


5 


0.313 


0.712 


0.081 


P29 


4 


0.250 


0.579 


4 


0.313 


0.591 


3 


0.500 


0.486 


4 


0.250 


0.659 


4 


0.813 


0.236 


0.046 


P30 


7 


0.500 


0.768 


5 


0.188 


0.623 


17 


0.125 


0.925 


6 


0.250 


0.750 


11 


0.375 


0.857 


0.000 


P31 


2 


0.938 


0.063 


2 


0.875 


0.121 


2 


0.938 


0.063 


1 


1.000 


0.000 


2 


0.938 


0.063 


0.066 


P35 


6 


0.375 


0.786 


7 


0.063 


0.702 


6 


0.313 


0.714 


7 


0.188 


0.756 


6 


0.375 


0.655 


0.794 


P41 


6 


0.375 


0.849 


7 


0.063 


0.726 


6 


0.250 


0.843 


6 


0.125 


0.827 


8 


0.250 


0.798 


0.000 


P44 


5 


0.000 


0.768 


4 


0.000 


0.599 


5 


0.125 


0.730 


5 


0.000 


0.726 


5 


0.438 


0.508 


0.000 


P45 


6 


0.133 


0.821 


6 


0.286 


0.725 


8 


0.063 


0.819 


5 


0.438 


0.736 


6 


0.273 


0.823 


0.000 


P47 


2 


0.875 


0.315 


2 


0.813 


0.272 


3 


0.688 


0.401 


3 


0.625 


0.325 


3 


0.750 


0.232 


0.316 


Average 


6 


0.389 


0.670 


5 


0.334 


0.557 


6.6 


0.359 


0.660 


5 


0.331 


0.616 


6 


0.374 


0.632 





Note: A = number of alleles; H e = expected heterozygosity; H Q = observed heterozygosity; HWE = Hardy- Weinberg equilibrium; N = sample size. 
a All values are based on 98 samples located in Xiaoxing'anling, Zhangguangcailing, Changbai Mountain, Laoyeling, and Wanda Mountain (Appendix 1). 
b Values shown are % 2 test results for each locus across all populations sampled. 
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Appendix 1. Geographic localities and sample sizes of the Juglans mandshurica populations in this study. Voucher specimens were deposited at the 
Northeast Forestry University Herbarium (Heilongjiang; NEFI). 



Locality 


Latitude (°N) 


Longitude (°E) 


Altitude (m) 


n 


Population code; herbarium voucher accession code 


Xiaoxing'anling 


46.8927 


128.4209 


275.0 


19 


XXAL; jmXXAL2010-Yu 


Zhangguangcailing 


44.0133 


127.0052 


357.5 


20 


ZGCL; jmZGCL2010-Yu 


Changbai Mountain 


41.1009 


126.1315 


332.0 


20 


CBS;jmCBS2010-Yu 


Laoyeling 


43.0405 


130.3383 


364.0 


19 


LYL; jmLYL2010-Yu 


Wanda Mountain 


44.5802 


131.1351 


115.0 


20 


WDS;jmLYL2010-Yu 



Note: n = sample size. 
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